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Abstract 

Purpose: To outline the evolution of computerized brachytherapy treatment planning in the United States through 
a review of technological developments and clinical practice refinements. 

Material and methods: A literature review was performed and interviews were conducted with six participants in 
the development of computerized treatment planning for brachytherapy. 

Results: Computerized brachytherapy treatment planning software was initially developed in the Physics De- 
partments of New York's Memorial Hospital (by Nelson, Meurk and Baiter), and Houston's M. D. Anderson Hospital 
(by Stovall and Shalek). These public-domain programs could be used by institutions with adequate computational 
resources; other clinics had access to them via Memorial's and Anderson's teletype-based computational services. 
Commercial brachytherapy treatment planning programs designed to run on smaller computers (Prowess, ROCS, 
MMS), were developed in the late 1980s and early 1990s. These systems brought interactive dosimetry into the clinic 
and surgical theatre. 

Conclusions: Brachytherapy treatment planning has evolved from systems of rigid implant rules to individualized 
pre- and intra-operative treatment plans, and post-operative dosimetric assessments. Brachytherapy dose distributions 
were initially calculated on public domain programs on large regionally located computers. With the progression of 
computer miniaturization and increase in processor speeds, proprietary software was commercially developed for 
microcomputers that offered increased functionality and integration with clinical practice. 
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Purpose 

The modern permanent prostate brachytherapy pro- 
cedure, introduced three decades ago [1], features 
transperineal deposition of radioactive sources under 
ultrasound-guidance. Superior disease control has been 
achieved due to accurate enactment of an implant plan 
that was developed either pre- or intra-operatively with 
computerized treatment planning system (TPS). While 
the general principles of brachytherapy treatment plan- 
ning have been covered elsewhere [2], the current paper 
reviews the evolution of computerized treatment plan- 
ning for brachytherapy in the United States. There were 
of course developments outside the U.S. during this 
same era, and the interested reader may learn more on 
these advances through examining the following refer- 
ences [3-6] . Although the authors recognize the contribu- 
tions of European and Asian workers to early advances 
in brachytherapy dosimetry, the scope of this paper is 
limited to the authors' area of our expertise - develop- 
ments in the U.S. 



Material and methods 

This study is based upon a literature review and origi- 
nal interviews, supported by contemporaneous literature, 
that were obtained from the following physicists and 
engineers: Lowell Anderson (March 2005), Stephen Baiter 
(June 2009 and February 2013), Philip Heintz (December 
2012), Ted Jackson (October 2012), Michael Saylor (Au- 
gust 2013), Marilyn Stovall (February 2013). 

Working in the dark 

In the earliest days, brachytherapists relied on atlas- 
es, tables, and experience to determine source distribu- 
tion, implant strength, and duration. By the middle of 
the 20 th century, interstitial brachytherapy (more proper- 
ly, br achy radium [7]) dosimetry relied on two source distri- 
bution systems, which were developed in the 1930s: Pat- 
erson-Parker and Quimby. The Paterson-Parker System, 
developed in Manchester's Holt Radium Institute, spec- 
ified an inhomogenous spatial distribution of sources to 
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achieve a relatively homogenous dose distribution [8,9]. 
The system developed by physicist Edith Quimby at New 
York's Memorial Hospital stipulated a homogenous dis- 
tribution of sources to generate an inhomogenous dose 
distribution [10,11]. Both systems utilized pre-planning 
(i.e., preceding clinical implantation) to determine the 
strength and geometric distribution of implanted sourc- 
es needed to deliver a specified dose (typically the mini- 
mum dose) to a target. 

The procedure used at New York's Memorial Hospi- 
tal during the 1950s and 1960s is described as an example. 
The brachytherapist was called to the operating theater to 
implant a tumor that had just been deemed unresectable 
(a common occurrence in an era without advanced im- 
aging). He had, at hand, a caliper and a supply of 222 Rn 
seeds. The total strength of implanted sources was calcu- 
lated by the following process: the three dimensions of 
the target were measured with a ruler and the volume 
calculated by using the appropriate formula for the tar- 
get shape. The maximally tolerated dose for the target 
volume was determined by consulting a Memorial chart 
(larger volumes were found to tolerate lower doses). 
Next, the Paterson-Parker or Quimby tables were used 
to determine the total milligram-hours of radium (i.e., 
mg-h 226 Ra, an obsolete unit of source strength) necessary 
to achieve the dose for that volume; the result was then 
multiplied by the appropriate Paterson-Parker elonga- 
tion-correction factor. Since most Memorial implants of 
that era used 222 Rn seeds, the mg-h 226 Ra value had to be 
divided by the 222 Rn mean life (i.e., 132 h) to determine 
the total implant strength. This value was divided by the 
current strength of each seed to determine the number of 
seeds to be implanted. Finally, a decision had to be made 
as to how to spatially distribute the seeds [12]. 

This process proved too cumbersome (especially as 
the patient's chest or abdomen had already been opened), 
so the following shortcut was implemented: the three 
measured dimensions were averaged. To determine the 
222 Rn source strength to implant, this average dimension 
was multiplied by 10 (an empirically derived factor, spe- 
cific to 222 Rn); a factor of 5 was used when 125 I seeds were 
introduced [12]. Thereafter, a Quimby-style (homoge- 
nous distribution of sources) implant was performed, 
which was appropriately fast and simple. Physicist Low- 
ell Anderson further simplified the process through the 
introduction of nomograms for 125 I and 103 Pd [13,14]. 

While the use of atlases, tables, and shortcuts sys- 
tematized the process, the implants were often serious- 
ly flawed. It was difficult (even for expert brachythera- 
pists) to reproduce the 'ideal' source geometry specified 
by the systems, and the achieved 'minimum' target dose 
typically fell well below the prescribed dose [15]. Larg- 
er tumors were frequently underdosed, a problem that 
was exacerbated by the introduction of 125 I. Initially, only 
the inverse-square law was used in calculations; this was 
not an issue for 226 Rn, 222 Rn, 198 Au, or 192 Ir due to their 
high-energy y-ray emissions, but was inadequate when 
calculating the attenuation of low-energy photons emit- 
ted by radionuclides such as 125 I [16]. Underdosing was 
not recognized because post-implant dosimetry was not 



initially performed. Manual calculation of dose at more 
than a few points was tedious, especially for a seed im- 
plant where the position of more than 50 sources had to 
be identified and their dose distributions had to be plot- 
ted and summed. 

What was needed was a system to identify a large 
number of sources, then calculate, sum, and display 
a spatial representation of the dose distribution over the 
entire volume of interest (not just the periphery). Com- 
puters have the capacity to calculate dose distribution 
rapidly enough, so that timely adjustments could be 
made post-operatively (for afterloaded implants), or even 
intra-operatively (for permanent implantation). This his- 
torical review explores the introduction of computer do- 
simetry to brachytherapy. 

Introduction of computerized dosimetry 
at Memorial Hospital 

Kia-Chi Tsien (1916-1993) introduced a TPS for mul- 
ti-field teletherapy using 'tabulating machines' in 1955 
[17,18] . The depth-dose characteristics, specified as a func- 
tion of field size and focus-to-skin distance, for all 60 Co, 
orthovoltage, and megavoltage equipment in use at Me- 
morial Hospital were stored in a library of punch cards 
(Fig. 1). The appropriate punch cards for the selected 
fields were entered into the tabulating machine, which 
summed their contributions and plotted the isodose lines 
relative to an external contour of the body. This system 
quadrupled the efficiency of dosimetrists, allowing selec- 
tion of the most favorable of several plans in the same 
time that it would have taken to run a single plan manu- 
ally. Richard Nelson and Mary Lou Meurk extended this 
work to perform the first computerized brachytherapy 
dose calculations in 1958 [19]. The positions of implanted 
sources were localized through stereo-shift radiographs 
and assigned a location at the nearest point on a three- 
dimensional Cartesian coordinate lattice with 5 mm spac- 
ing (Fig. 2). The location of each source was transferred 
to a punch card, and the tabulating machine summed 
the contribution of all sources to plot the dose distribu- 
tion. Output was represented as a matrix of points with 
1 cm spacing; isodose lines were drawn by hand. By 1961, 
the system was programed in the FORTRAN code for 
running on a time-shared model 7090 IBM (International 
Business Machines, Corp., Armonk, NY, USA) comput- 
er system, which was "sufficiently fast that its results 
have been used to modify interstitial implants before and 
during treatment" [20]. While the system was feasible 
for 226 Ra needles and 192 Ir ribbons in afterloading tubes, 
localization of individual seeds was "possible, but not 
always practicable" as the program accommodated an 
"equally spaced array of seeds" [20]. 

Stephen Baiter (born in 1940) revised the system 
shortly after joining the Memorial Physics staff in 1963. 
He had learned how to program the UNIVAC (UNIVer- 
sal Automatic Computer) by Eckert-Mauchly Computer, 
Corp. (Philadelphia, PA, USA) while he was an under- 
graduate student at the Polytechnic Institute of Brook- 
lyn (B.S. in Physics, 1961), and had access to Colum- 
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Fig. 1. A paper punch card used for computerized dose calculations at Memoria 
obtained from the British Institute of Radiology 
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Hospital from Tsien 1958 [14]. Permissions 



bia University's computer laboratory during graduate 
work (M.S. in Radiologic Physics, 1963). He learned the 
FORTRAN II computer programming language from 
a manual while commuting to work by subway. It was 
his intention to write a program for Memorial's model 
CDC-160A (Control Data Corp., Minneapolis, MN, USA) 
4k 12-bit minicomputer to spatially display the dose 
distribution in any plane for a large seed implant, to 
calculate average and minimum doses, and to generate 
dose-volume histograms [15,21]. The program had to be 
broken into modules that ran sequentially because com- 
puter memory could not store the entire program. First, 
seed location was determined by the stereo-shift method. 
This was a tedious task because implants often contained 
more than 50 seeds (Fig. 3). After the source locations 
were entered, dose was calculated to points on a matrix. 
Dose calculation points in the vicinity of the target vol- 
ume were typically 1 cm apart, so there were 1,000 dose 
points for a 10 cm 3 volume. It took about 0.1 s to calculate 
dose to a point from a single source, 100 s to calculate 
the dose to all dose points from a single source, and over 
an hour to calculate the sum of contributions from all 
sources to all dose points. A printout of dose to points 
on any plane was generated, but isodose lines were still 
drawn by hand (Fig. 4). The situation improved when the 
Memorial Physics Department installed a stereographic 
reference frame (Fig. 5) to maintain fixed geometry for 
stereoshift or orthogonal radiographs, a digitizer to enter 
seed location, a program to automatically digitize seed 
location, and a model 1800 IBM 16k computer that in- 
creased calculation speeds. 

Although the 'Memorial Implant Dosimetry Appli- 
cation System' ('MIDAS') was freely shared, few institu- 
tions had the computing power to utilize it. Beginning in 
1967, the Memorial Physics Department offered a com- 
putational service, whereby other institutions submitted 
implant source strength and spatial distributions. Dose 




100% 



Fig. 2. A) Dose summation of four brachytherapy sources 
onto a plane, and B) creation of isodose distributions from 
Nelson and Meurk 1958 [15]. Permissions obtained from 
the Radiological Society of North America 

distributions were calculated using a model 360/67 IBM 
central processor on the General Electric (GE, Corp., 
Fairfield, CT, USA) commercial timeshare service, and 
results were returned to the participating institution by 
teletype [22]. Initially, six institutions utilized the ser- 
vice [23], but more than 60 hospitals participated by 
1972 [24]. 
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Fig. 3. Anterior-posterior drawing from 1955 of the posi- 
tion of 222 Rn seeds within the right lung for brachytherapy 
dose calculation using the MIDAS TPS. This image was 
provided courtesy of Stephen Baiter 



Fig. 4. Printout of 222 Rn point doses on a plane by the 
MIDAS TPS, with isodoses delineated by hand. This im- 
age was provided courtesy of Stephen Baiter 



Introduction of computer dosimetry 
at M. D. Anderson Hospital 

Similar events transpired concurrently in Houston. 
Unlike at New York's Memorial Hospital where 222 Rn 
seeds were commonly implanted, the radiotherapists at 
M. D. Anderson Hospital preferred 226 Ra needles and 
198 Au seeds. Upon graduation from Baylor University 
in 1951 with a B.S. degree in Mathematics and Zoology, 
Marilyn Stovall (born in 1931) answered a newspaper ad- 
vertisement for employment at Anderson. Robert Shalek 
(1923-2007), acting Head of Physics, offered her a dosim- 
etrist position with on-the-job training if she would corn- 




Fig. 5. The stereographic reference frame developed by 
the Memorial Physics Department to maintain a fixed 
geometry for stereoshift or orthogonal radiographs. This 
image was provided courtesy of Stephen Baiter 



mit to remaining there for at least a year. She accepted his 
offer and has been there for more than six decades. 

Her initial responsibilities included machine calibra- 
tion, treatment planning, and checking dose delivery us- 
ing a Sievert chamber in the radiation field. As custodian 
of the 226 Ra sources, she was present at the time of im- 
plantation and also when the sources were removed (i.e., 
at any hour of the day or night). After several years, her 
position was advanced to 'Applied Mathematician', and 
she was asked to write computer code for brachytherapy 
dose calculations. She learned the FORTRAN computer 
programming language at the University of Houston's 
night school, and wrote a program with the assistance of 
a professor at Texas A&M; punch cards were transported 
between Houston and College Station, TX on passenger 
bus. The program was ready in 1960, but the only usable 
computer, a RAMAC-model IBM, was in the business of- 
fice. The program was first used to calculate post-implant 
dosimetry for radium needle implants using orthogonal 
x-ray images, then still commonly referred to as 'roent- 
genograms'. Source localization was done by transverse 
tomography. The computer determined the square of the 
distance to points on a 1 mm grid, assigned doses based 
on a tabular lookup, and summed the doses [25]; isodose 
curves were drawn manually. It took an hour for each 
plane to be calculated, at which point a new set of cards 
would be fed into the computer. As dose points were 
spaced at 0.5 cm (there were 8 times as many dose cal- 
culation points per volume as had been utilized at Me- 
morial), the entire process initially took 10 hours. As the 
computer was used for billing during the day, Stovall set 
up a cot in the business office and ran post-implant do- 
simetry at night. This situation continued until the An- 
derson Physics Department acquired its own faster corn- 
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puter, a model 1620 IBM, which could perform the same 
task within minutes. 

In a retrospective analysis of implants, Stovall had de- 
termined that most recurrences and complications (e.g., 
tissue necrosis) occurred in regions of under- or over- 
dosage [26], respectively. Timely dose calculation, with 
rapid depiction by plotter or cathode ray tube, allowed 
compensation for less than 'ideal' geometry by varying 
the timing of needle removal for interstitial implants and 
through source selection for afterloading implants [27]. 

Similar to Memorial's MIDAS dosimetry program, 
Anderson's RADCOMP (RADiation COMPutation) do- 
simetry program, was in the public domain [28]. As few 
institutions had the computing capacity to utilize it, the 
Anderson Department of Physics also had a thriving 
multi-institutional computational service. 

Advances in software algorithms 

After the introduction of computerized brachythera- 
py treatment planning, advancements beyond the point- 
source approximation were introduced to account for 
the elongated physical distribution of the radionuclide 
[29,30]. Cooperating with Shalek at the Anderson Depart- 
ment of Physics, Meisberger et al. addressed the issue of 
radiation attenuation and scatter in the medium through 
dose correction factors that were expressed as simple 
polynomial functions [31]. While initially used only for 
high-energy photon-emitting sources such as 198 Au, 192 Ir, 
137 Cs, 226 Ra, and 60 Co in the late 1960s, this work was later 
extended to low-energy photon-emitting sources such as 
125 I, where the radiological influence of medium is more 
pronounced [32]. Another important influence on clinical 
dose calculations in this era was Monte Carlo methods 
for radiation transport simulations, where the first full 3D 
model of a brachytherapy source was performed in 1971 
[33]. Results from these types of Monte Carlo simulations 
were simplified for integration into existing planning al- 
gorithms and development of future systems. 

Development of proprietary software 

Peter Wootton (1924-2004), departed Anderson to join 
the Physics faculty at the University of Washington. One 
of his doctoral students, Philip Heintz (born in 1943), had 
learned computer programming as an undergraduate en- 
gineering student at California Polytechnic State Univer- 
sity from 1961 to 1965. Together with physicist Douglas 
Jones, they used elements of the Memorial and Anderson 
public domain programs to generate the TPS at the Uni- 
versity of Washington. Heintz took the program with him 
when he went into private practice in California in 1971. 

Heintz joined Ralph Worsnop in 1985 to develop MARx 
Plan, a program written in the PASCAL computer pro- 
gramming language to be run on a personal computer (PC). 
The partnership with Worsnop dissolved in 1987 and Heintz 
had the program re-written in the Microsoft C (Microsoft, 
Corp., Redmond, WA, USA) computer programming lan- 
guage to be run on the IBM AT platform. His new company 
SSGI (Small Systems Group, Inc., Chico, CA, USA) intro- 



duced the Prowess 2000 TPS software in 1988. In addition 
to teletherapy planning, it had modules for brachytherapy 
with ^Ra and 137 Cs tubes, 192 Ir and 125 I seeds, and 182 Ta 
wire. In the early 1990s, the program was rewritten for the 
Windows (Microsoft, Corp., Redmond, WA, USA) operat- 
ing system and was renamed Prowess 3000. 

Another commercial brachytherapy TPS was the 
ROCS program (Radiation Oncology Computer Systems, 
Inc., Carlsbad, CA, USA) introduced in 1988 by physicist, 
founder, and president David Neblett. It was written in the 
BASIC computer programming language for the PC and 
included an ultrasound-guided prostate brachytherapy 
package. 

William Saylor was trained as a medical physicist at 
the Naval Postgraduate School in Monterey, CA in 1962. 
In 1974, he introduced Therpac (Pinecliff Associates, 
Wilmington, NC, USA), the first complete TPS written 
for the model HP 9830 A (Hewlett-Packard, Corp., Palo 
Alto, CA, USA) microcomputer. His son, Michael Saylor, 
joined him in 1981 and rewrote the program in Microsoft 
C computer programming language for the PC; Ther- 
pac-PC was introduced in 1986. In 1993, a transperineal 
ultrasound-guided implant (TUI) module (Therpac Plus 
6.5 B3Dtui) was introduced. Within a few years, the mod- 
ule could incorporate ultrasound images for pre-planning 
and CT images for post-implant dosimetry. William Say- 
lor retired in 1996, and sold the company to Multimedia 
Medical Systems (MMS), Inc. (Charlottesville, VA, USA), 
a telemedicine company specializing in video-conferenc- 
ing and image-sharing. Michael Saylor rewrote the TUI 
module with graphical-user interface for the Windows 
operating system, introduced in 1997. MMS's chief engi- 
neer, Ted Jackson (Biomedical Engineering M.Sc. in 1990 
and Ph.D. in 1994 from the University of Virginia) collab- 
orated with the physicians and physicists at the Mount 
Sinai Hospital in New York to develop a user-defined 
dose-optimization module, and the facility for interactive 
intra-operative planning with real-time display of an im- 
planted dose distribution. Two years later, the MMS pro- 
gram was acquired by Varian Medical Systems, Inc. (Palo 
Alto, CA, USA) and renamed 'VariSeed™'; it has been 
the dominant program in use in the U.S. for interstitial 
brachytherapy treatment planning for over a decade. 

Conclusions 

Computerized brachytherapy treatment planning 
software was initially developed in the Physics Depart- 
ments of Memorial Hospital (by Nelson, Meurk, and 
Baiter) and M. D. Anderson Hospital (by Stovall and Shal- 
ek). These systems were public-domain, and available 
to smaller hospitals through multi-institutional compu- 
tational services. Commercial brachytherapy treatment 
planning programs, designed to run on microcomputers 
computers, were introduced in the late 1980s and early 
1990s. These systems brought computer brachytherapy 
dosimetry into the clinic and surgical theatre. 
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